Pomeron structure functions from HERA to Tevatron and LHC 
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The proton diffractive structure function F^'^' measured in the HI and ZEUS experiments at 
HERA are analysed in terms of perturbative QCD in the perspective of the QCD extrapolation 
to the Tevatron and the LHC. It is shown that both data sets can be well described by a QCD 
analysis in which point-like parton distributions, evolving according to the next-leading DGLAP 
equations, are assigned to the leading and sub-leading Regge exchanges. For present data from HI 
and ZEUS the gluon distributions are found to be quite different and we give the corresponding sets 
of quark and gluon parton distributions for the Pomeron, extracted from the two experiments. An 
extrapolation to the Tevatron range is compared with CDF data on single diffraction. Conclusions 
on factorization breaking between HERA and Tevatron critically depend on whether HI (strong 
violation) or ZEUS (compatibility at low f3) fits are taken into account. Using the double Pomeron 
formulation in central diffractive dijet production we show that the Tevatron mass fraction is much 
^~>) ■ sensitive to the high (3 tail of the gluon in the Pomeron, suggesting a new way of handling the 

' otherwise badly known gluon distribution in the Pomeron. Extrapolation of the fits to very high 

are given since they will be relevant for QCD and diffraction studies at the LHC. 

. I. INTRODUCTION 

Since years, the Pomeron remains a subject of many interrogations. Indeed, defined as the virtual colourless carrier 
of strong interactions, the nature of the Pomeron is still a real challenge. While in the perturbative regime of QCD it 



o 



> 



can be defined as a compound system of two perturbatively Reggeized gluons [1] in the approximation of resumming 
^ the leading logs in energy, its non-perturbative structure is basically unknown, 
psj ■ However, in the recent years, an interesting experimental investigation on "hard" diffractive processes led to a new 
CN| ' insight into Pomeron problems. It is now experimentally well established at HERA [2-4] that a substantial fraction 
of ep events is contributablc to diffraction, i.e. color singlet exchange, initiated by a highly virtual photon. Starting 
■ with the pioneering theoretical work of Ref. [5] , the idea of a point-like structure of the Pomeron exchange opens the 
JS^ I way to the determination of its parton (quark and gluon) distributions, where the Pomeron point-like structure can 
-««^ ■ be treated in a similar way as (and compared to) the proton one. Indeed, leading twist contributions to the proton 
' diffractive structure functions can be defined by factorization properties [6] in much the same way as for the full 
^p-i. proton structure functions themselves. As such, they should obey DGLAP evolution equations [7], and thus allow for 
Q_i' perturbative predictions of their evolution. 

D , In Ref. [8], we have shown that the data are well described by a QCD analysisin which point-like parton distributions, 
' evolving according to the DGLAP equations, are assigned to the leading and sub-leading Regge exchanges. The gluon 
J> I distributions were found to be quite different for HI and ZEUS. Then, in Ref. [8], we have derived sets of quark and 
gluon parton distributions for the Pomeron, and predictions for the charm and the longitudinal diffractive structure 
function from the QCD fit. An extrapolation to the Tevatron range was also compared with CDF data on single 
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^ 1 diffraction. In Ref. [9], we have asked the question whether the quark and gluon distributions in the Pomeron 
■ " ' obtained from QCD fits to hard diffraction processes at HERA could be dynamically generated by QCD evolution 
from a state made of valence-like gluons and sea quarks as an input. By a method combining backward Q^-evolution 
for data exploration and forward Q^-evolution for a best fit determination, we have found that the diffractive structure 
functions published by the HI collaboration at HERA could be described by a simple valence-like input at an initial 
low scale. The same property was not achieved when using ZEUS data. 

These previous analyses were based on diffractive deep inelastic scattering (DIS) cross-sections determined at HERA 
[2,3]. It can be seen as a first stage in studying aspects of the QCD fit technique applied to diffraction, as for example 
the flexibility of the parton distributions parameterisations. 

In the present paper, we produce new sets of diffractive parton distribution functions (DPDFs) extracted from the 
presently published data sets, which consist in the already mentionned set for HI [2] and the latest results from ZEUS 
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difFractive DIS cross-sections [4]. The results on DPDFs arc close to those of Ref. [8]. We give for the first time a 
new method to evaluate the uncertainty on the gluon density at large values of /?, which is a key point for any further 
discussion of the influence of DPDFs at Tevatron or LHC. 

Then, we present the main result of this article, which is to show how these diffractive parton densities enter in 
the simulations of the dijet mass fraction at the Tevatron in central difFractive dijet production within the double 
Pomeron formulation. We show that this quantity is much sensitive to the high (3 tail of the gluon in the Pomeron. 
Our analysis allows to provide the extrapolation of the fits to very high since they will be relevant for QCD and 
hard diffraction studies at the LHC. 



II. EXTRACTION OF PARTON DISTRIBUTIONS IN THE POMERON AND QCD FITS 

A. Formulation 

It is well known [2,8] that the diffractive structure function F^^^\ measured from DIS events with large rapidity 
gaps, can be investigated in the framework of Regge factorization and expressed as a sum of two factorized contribu- 
tions corresponding to a Pomeron and secondary Reggeon trajectories : 

F^^'\q\P,xp) = /p/p(a;p) if (Q^/3) + /r/p(xp) F^{Q\f3) . (1) 

In this parameterisation, can be interpreted as the Pomeron structure function and F^ as an effective Reggeon 
structure function, with the restriction that it takes together into account various secondary Regge contributions which 
can hardly be separated. The Pomeron, /p/p, and Reggeon, fu/p, fluxes are assumed to follow a Regge behaviour 
with linear trajectories ap^u{t) = <^ip,u{0) + <^'p^u^, such that 

/p/p,R/p(a;p) = / dt (2) 

where |trnm| is the minimal kinematically-allowed value of |t| and tcut = ^1 GeV^ is the limit of the measurement. 
The values of the t-slope parameters are taken from hadron-hadron data (ap = 0.26 GeV^^, = 0.90 GeV~^, 
Bp = 4.6 GeV^^, Br = 2.0 GeV~^), since the data are not precise enough to determine them [2]. 

We assign parton distribution functions to the Pomeron and to the Reggeon, considered as coumpound states 
of fundamental quarks and gluons. A simple and well-known prescription of QCD factorisation is that the parton 
distributions of both the Pomeron and the Reggeon are parameterized in terms of non-perturbative input distributions 
at some low scale which was adopted here to be Qq = 3 GeV^. The pion structure function [13] is assumed to be also 
valid for the sub-leading Reggeon trajectory with a free global normalization to be determined by the data ^ . 

For the Pomeron, a quark flavour singlet distribution, zS{z, Q"^) = u + u + d + d + s + s, and a gluon distribution, 
zG{z,Q'^), are parameterized in termsof coefficients C^^^ and cj*^^ at Qq = 3 GeV^ as it was done in Ref. [2] such 
that 



zS{z, = Ql) = Cf . Pj{2z - 1) j . exp i^-^j (3) 
.G(z,Q^ = Q2)= j^cf).p^.(2^_l)j -exp^-^) (4) 



where z = Xi/p is the fractional momentum of the Pomeron carried by the struck parton, Pj{C) is the j*'* member in 
a set of Chebyshev polynomials, chosen such that Pi = 1, P2 = C) and Pj+i{() = 2(Pj{() — Pj-i{(). In the following. 



^We checked that changing the pion structure function by some amount (20%) does not change significantly the parton 
distributions. For even better checking, another pion structure function [14] has also been used without any sizable modification 
of the results. 
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we present the results in terms of the overall quark density with the hypothesis that u = u = d = d = s = s, which 
means that each light quark density is equal to zS/6. 

A sum of n = 3 orthonormal polynomials is used so that the input distributions arc free to adopt a large range of 
forms for a given number of parameters. Any bias towards a particular solution due to the choice of the functional 
form of the input distribution is therefore minimized. 

In Ref. [8], we have used the following fixed Pomeron intercepts for HI and ZEUS data sets, Q:p(0) = 1.20 ± 0.09 
and ap(0) = 1.13 ± 0.04 respectively, as well as qr(0) = 0.62 ± 0.03 for the Reggeon intercept in case of HI (there is 
no need for Reggeons for ZEUS data). In the following, we introduce the Pomeron intercept as a free parameter in 
the QCD fit. 

B. Data sets 

Only data points with virtuality > 3 GeV^, diffractive mass at the photon vertex Mx > 2 GeV and rapidity 
y < 0.45 arc included in the fit in order to avoid large higher-twist effects and the region that may be most strongly 
affected by a non-zero value of the longitudinal over transverse ratio R. 

The fits include 179 data points for HI data [2] with 8 parameters (3 for the sea quark density, 3 for the gluon 
density and 1 for the normalization of the Reggeon contribution and 1 for the Pomeron intercept). Thus, we have 171 
degrees of freedom with a of 203.2 using statistical errors only. In Ref. [8], only 161 data points were included in 
the QCD fit procedure since an additional cut for /3 < 0.65 was considered in the analysis. In the new version of the 
fits, we gain in sensitivity with a larger number of data points, in particular on the gluon density. For example, the 
third parameter for the gluon density C3 was not constrained in Ref. [8] with a value of 0.01 ± 0.04 and it is found 
to be zero in the new analysis. Then, for HI data, 3 parameters for quarks and 2 parameters for the gluon densties 
are enough to give a proper description of the cross section measurements through the QCD fit procedure. Also, 
new results and parameters obtained in Ref. [8] are in perfect agreement for the quarks within the quoted statistical 
uncertainties. 

For ZEUS data [4], the fits include 102 data points with 7 parameters (the parameter corresponding to the normal- 
ization of the Reggeon contribution is absent). We have 95 degrees of freedom with a of 118.4 using statistical errors 
only. Note that ZEUS measurements arc realised for diffractive mass at the proton vertex My < 2.3 GeV and squared 
transfer- momentum t < 1 GeV^, whereas HI gives diffractive cross sections for My < 1.6 GeV and t < 1 GeV^. We 
have converted ZEUS data points to the same My range by multiplying ZEUS values by the factor 0.77 [4]. This 
factor takes into account the conversion from My < 2.3 GeV to My < 1.6 GeV for the ZEUS measurements, which 
reflects mainly the fact that the proton-dissociation is reduced for the range of HI data set. The uncertainty on this 
number is of course large (of about ±0.1 [4]), but we do not take it into account in the following as we do not propose 
a common fit of HI and ZEUS data sets, which would require a relative normalisation error. 

On Fig. 2, structure function measurements from both data sets are shown. In this last figure, ZEUS data points 
have been multiplied by the factor 0.77 as explained above, and extrapolated to HI bin centers. This extrapolation 
is done just for the plot, ZEUS data points have been included in the QCD fits at their exact kinematic values. 

To get solvable evolution equations, the parton distribution functions must approach zero as 2; ^ 1. This is achieved 
by introducing in Eqs. (4) the exponential term with a positive value of the parameter a. Unless otherwise indicated, 
in the following fits a is set to 0.01 such that this term only infiuences the parametrisation in the region z > 0.9. 
This term is only present to ensure the convergence and plays a role in a domain where we do not include data (for 
instance, at /3 =0.65, its value is equal to 0.97). 

The functions zS and zG are evolved to higher using the next-to-leading order DGLAP evolution equations 
(with the code of Ref. [18]), with as(M|) = 0.118 (in Ref. [8], a lower value of as(M|) has been considered). The 
contribution to i^2^(/3, Q^) from charm quarks is calculated in the fixed flavour scheme using the photon-gluon fusion 
prescription given in Ref. [19] as implemented in Ref. [20]. The contribution from heavier quarks is neglected. 

No momentum sum rule has been imposed because of the theoretical uncertainty in specifying the normalization of 
the Pomeron or Reggeon fluxes. It is also not clear that such a sum rule is appropriate for the parton distributions of 
a virtually exchanged state. After fitting, however, we observed that the sum rules were fulfilled within 10% despite 
being not required in the fits. 

C. Results 

Already from Fig. 2 we can get a good idea of the relative gluon densities between HI and ZEUS data sets. Indeed, 
we notice a good agreement at Q'^ = 7.5 GeV^ (and /? < 0.65) whereas both sets of measurements show discrepancies 
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at larger = 28. GeV^. ZEUS cross section values are below HI for this value, which means that scaling 
violations in case of ZEUS are lower. Thus, from this basic observation, we expect a lower gluon density in case of 
ZEUS data compared to HI. 

The resulting parton densities of the Pomeron extrapolated from the QCD fit procedure are presented in Fig.l 
(extrapolated till 10000 GeV^) and the parameter values are given in table I. 
A few comments are in order : 

• We notice that the quark contribution is smaller compared to the gluon one, and much smaller for HI prediction, 
confirming the standard results exposed for example in Ref. [8]. On Fig.l, the uncertainties of the parton 
distributions are not shown but they can be derived from table I, giving typically a 25% error for the gluon 
density. 

• In the range of the measurements, essentially below = 75 GeV^ for HI [2] and below = 55 GeV^ for 
ZEUS [4], the agreement between HI and ZEUS for quark densities is reasonnable. 

• Gluon densities for both data sets is quite different, especially at large /?. This last point has been addressed 
above from the simple observation of scaling violations and is more extensively discussed below. 

The result of the fit is presented in Fig. 2 together with the experimental values. The fit is plotted only in the 

(3 — Q"^ bins considered for the QCD analysis. In this last figure, as mentioned above, ZEUS data points have been 
converted to the same My < 1.6 GeV domain as HI (multiplying ZEUS measurements by a multiplicative factor 0.77 
as explained above) and extrapolated to HI bin centers. However, the QCD-fit procedure and extraction of parton 
distributions have been done exactly on ZEUS published data points [4]. The extrapolation to HI bin centers is 
applied only for the presentation of Fig. 2 and we have checked that these conversion factors are small (of the order 
a few per cent). We sec on this figure the good agreement of the QCD fit with the data points, which supports 
the validity of the description of the Pomeron in terms of partons following perturbative QCD dynamics, even if the 
accuracy of the gluon determination is still quite poor. This provides a very nice perspective for diffractive analyses 
of more recent data to be released, where better accuracy will lead to a more precise QCD analysis. 



parameters 


HI 


ZEUS 


ap 


1.19 ± 0.02 


1.13 ± 0.02 




0.21 ± 0.05 


0.38 ±0.02 


MS) 

^2 


0.04 ± 0.02 


-0.03 ± 0.01 




-0.14 ± 0.02 


-0.11 ± 0.01 




0.59 ± 0.40 


0.39 ± 0.10 




-0.01 ± 0.01 


-0.36 ± 0.05 




0.00 


0.05 ±0.02 


Nik 


14.11 ± 1.13 


0.00 



Table I- Pomeron quark and gluon densities parameters. Note that ' is found to be negligeable for HI data. 
Also, as mentionned above, there is no Reggeon contribution for ZEUS data. For HI data, the parameters are close 
to those of Ref. [8] within the quoted statistical uncertainties, even if we expect some differences as it is explained in 
the text : for example, more data points are considered and some small changes are done in the QCD fit procedure. 

D. Large z behaviour 

In order to analyze in more detail the large z behaviour of the gluon distribution zG{z, = Qq) and give a rough 
estimate of the systematic error related to our parametrizations, we consider the possibility to change the ansatz (4) 
by a multiplicative factor (1 — z)" . We shall consider i' to vary in the interval ±1 in order to allow for the still large 
indetermination of the gluon distribution. If we include this multiplicative factor (1 — z)" in the QCD fit analysis 
with statistical errors, we derive a value oi i' — 0.0 ± 0.6. Then, the variation in the interval ±1 considered here takes 
also into account a systematic effect of the order of the statistical uncertainty, which is correct. As it is explained 
more extensively below, the understanding of the large z behaviour is of essential interest for any predictions at the 
Tevatron or LHC in central dijets production. In particular, a proper determination of the uncertainty in this domain 
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of momentum is necessary and the method we propose in this article is a first estimate, that can be propagated to 

other measurements. 

Following Rcf. [30], one knows that starting with a given large z input function of the form (1 — zY° , the effect of 
leading order DGLAP evolution is to change the end-point exponent by 

(X-zY- - (\-zT-^ ^o^'^s{Q^) ^ (5) 

where n/ is the number of active flavors, up to a normalization factor. 

Hence, let us consider a given solution of the DGLAP evolution equation, such as our structiu'c function zG{z^ Q^), 
and its large z tail at Qq parametrized by (1 — zY°. Then, changing the high z behaviour of the input function (4) 
by (1 — z)^" will change the exponent y{Q^) by the same shift, namely 

y{Q^) -^yo- log a«(g') ±y = y{Q^) ± u . 

66 — ZUf 

Note that the same property is valid for the singlet quark distribution [30] . 

To be fully consistent with the NLO analysis, wo estimated the effect of next-leading DGLAP evolution on the 
modifications of the high z behaviour due to next-leading effects [31]. We found that the next-leading corrections to 
the overall high z exponent l/v{Q'^) is of the order of its inverse l/v{Q'^) < 20 % when v{Q'^) varies in the kinematical 
interval in consideration. This correction is negligeable in the conditions of our fits. 
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FIG. 1. Quark and Gluon distributions in the Pomeron.{}ll : full line, ZEUS : dotted line) as a function of z, the fractional 
momentum of the Pomeron carried by the struck parton, from the fit on HI and ZEUS data points with > 3 GeV'^. The 
parton densities are normalised to represent xp times the true parton densities multiplied by the flux factor at xp — 0.003. 
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FIG. 2. Comparison of data with the results of the QCD fits. Data points: HI (dashed line); ZEUS (full line). QCD fits: 
HI (crosses); ZEUS (diamonds). The ZEUS data points have been converted to the same My range as HI measurements and 
moved to the HI bin centers (see text). 
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FIG. 3. Gluon distributions of the Pomeron from the fit on HI data points alone with > 3 GeV^ . The parton densities 
are normalised to representor times the true parton densities multiplied by the flux factor atxp = 0.003. A change by 
v = ±.5 and ± 1 have been performed (see text). 
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FIG. 4. Gluon distributions of the Pomeron from the fit on ZEUS data points alone with > 3 GeV' . Same conventions 
as previous figure. 
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III. EXTRAPOLATION TO TEVATRON AND COMPARISON WITH CDF DATA 



The QCD fits we obtained from HERA data allow us to make direct comparisons for measurements at the Tevatron. 
It is quite interesting to be able to test directly the factorization breaking between HERA and the Tevatron using the 
measurements performed at both accelerators. We thus compare the extrapolations of the HI and ZEUS QCD fits 
to the recent CDF single diffractive jet cross-section measurement [11]. The result is given in Fig. 5. We note a large 
discrepancy both in shape and normalization between HI predictions and CDF data, clearly showing factorization 
breaking. However, the ZEUS fits are more compatible in normalization with the CDF measurement even if the shape 
is not described properly. We know that the gluon density from ZEUS is between 2 and 3 times smaller than the one 
from HI. The predictions for single diffractive production of jets at the Tevatron are thus expected to be different 
by an sizeable factor 2 to 3 between HI and ZEUS since they correspond to single Pomeron exchange. If the large 
statistical and systcmaticimccrtaintics on the gluon density are taken into account (about 50%for ZEUS, 25% for HI), 
ZEUS data are compatible with factorization at low /?. 

One should, however, question whether one has the right to extrapolate ZEUS results without introducing the 
additional Rcggcon component. Namely the CDF measurement is in a region in xp where the Reggcon contribution 
is important, in contrary to the ZEUS measurement. The error bar on the extrapolation from the fit for ZEUS is thus 
enhanced due to this uncertainty. Concerning the extrapolation of the HI fit and the comparison with CDF data, 
one notices that the CDF data lie primilarly at low values of /? where there are few HI data points. Moreover, in this 
kinematical domain, the HI data have the tendency to lie more in the high xp region where the Reggeon contribution 
is important and not well constrained by the fit. Thus, the extrapolation to the CDF domain suffers from large 
uncertainties. A combined fit using CDF data to constrain the low (3 region and the HERA data to constrain the 
high (3 domain would thus be of great interest. 

More precise data from HERA and detailed comparisons with the Tevatron are thus needed to study precisely 
factorization breaking between both experiments. The discussion of eventual higher twist contributions is clearly also 
valuable in order to reach conclusions on the shape at larger /3. It is thus important to get an accurate measurement 
of the gluon density in the Pomeron at low values of /? from HERA. Furthermore the Forward Proton Detector [29] 
installed by the DO collaboration for Run II will be of great help to get a direct measurement of the diffractive 
structure functions. 



IV. POSSIBLE MEASUREMENTS AT THE TEVATRON 

The CDF and D0 experiments at the Tevatron can measure directly the dijet mass fraction, defined as the ratio 
of the dijet mass to the total diffractive mass measured either in roman pot detectors or in the main D0 or CDF 
detectors. 

Let us first notice that we normalised the diffractive dijet cross section using directly the CDF run I measurement 
previously which takes into account the factorisation breaking between the Tevatron and HERA [8,11]. It has already 
been mentionned that this measurement is quite sensitive to the Pomeron structure in quark and gluons [8]. We 
display in Fig. 6 the dijet mass fraction for the two different gluon distributions described previously from HI and 
ZEUS data after requiring two jets in the central detector with a transverse momentum greater than 25 GeV. We 
indeed see that the difference between the gluon distributions from HI and ZEUS at HERA can be probed using 
Tevatron data, provided the survival gap probability (correction which is due to soft interactions) is constant over 
the full kinematical range. To better show the possibility of such a measurement, we give in Fig. 7 the dijet mass 
fraction when the p is tagged in the dipole roman pot detector from the D0 collaboration. The dipole detectors show 
a good acceptance for t close to 0. It would also be possible to require double tagged events (the antiproton in the 
dipole roman pot detector and the proton in the quadrupole one), but the acceptance of the quadrupole detectors 
only extends to i ~ 0.8 GeV^ which cuts a lot of events. 

On the other hand, it is also possible to probe the high-/? gluon density in the Pomeron at the Tevatron again 
through the measurement of the dijet mass fraction (or the total diffractive mass). In Fig8, we givethe dijet mass 
fraction using the shape of gluon distribution in the Pomeron coming from HI. The sensitivity to the uncertainty on 
the gluon distribution at high P is indicated on that figure by multiplying the gluon distribution by (1 — /?)'', which 
enhances or decreases the high /3 gluon distribution. It is quite important to be able to constrain this distribution 
since it is a direct background to an eventual exclusive signal at high /3. We notice that it is indeed possible to 
constrain the high-/? gluon using Tevatron data. At the LHC, it will be possible to constrain better the high-/? gluon 
density using higher mass objects, for instance in tt inclusive diffractive production. 
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FIG. 5. P distributions of CDF data compared with extrapolations from the diffractive parton densities. The partons densities 
are extracted from F^^^^ measurements by HI and by ZEUS collaborations 
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FIG. 6. Results for Dijet mass fraction at generator level. The HI and ZEUS parton distributions in the Pomeron serve as 
inputs (see text). 
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FIG. 7. Results for Dijet mass fraction in the D0 dipole roman pot acceptance at generator level. One considers tagging a p 
(see text). Same convention as previous figure. Note that the CDF results would be similar. 




dIjet mass fraction 

FIG. 8. Results for Dijet mass fraction for the HI parton distributions as a function of the high-B tail. The high-,5 component 
of the gluon density in the Pomeron is modified ed by a factor (1 — Z?)*" to show the dependence of the dijet mass fraction on 
the parameter v. (see text). 
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V. CONCLUSION 



We have shown that the proton diffractive structure function F2 measured in the HI and ZEUS experiments 
at HERA can be well described by a pcrturbativc QCD analysisin which fundamental quark and gluon distributions, 
evolving according to the next-leading DGLAP equations, are assigned to the leading and sub-leading Regge exchanges. 

The gluon distributions have been found to be quite different when extracted from present published data for HI 
and ZEUS. An extrapolation to the Tevatron range have been compared with CDF data on single diffraction, leading at 
this stage to quite different conclusions on factorization breaking between HERA and Tevatron depending on whether 
HI (strong violation) or ZEUS (compatibility at low /?) fits are used. The discrepancies of the gluon densities between 
both data sets have also been explained when looking simply at the scaling violations, the QCD fit procedure giving a 
quantitative estimate of these differences. We have presented a new method to evaluate the uncertainty on the gluon 
density at large values of /?, which is a key point for any further discussion of the influence of DPDFs at Tevatron or 
LHC. 

In double diffractive dijet production in the central rapidity region, within the double Pomeron formulation, we have 
found that the Tevatron mass fraction is very sensitive to the high [3 tail of the gluon in the Pomeron. Extrapolation 
of the fits to very high Q-^ have been given since they will be relevant for QCD studies at the LHC. 

Hence, we have shown that it is possible to constrain the high-/3 gluon using Tevatron data, which may be of decisive 
interest for any future measurement in central diffractive dijet topologies. 

At the LHC, it will be possible to constrain better the high-/? gluon density using higher mass objects, for instance 
in ti inclusive diffractive production. 

We hope that our parametrisations, and the ones easily obtainable exactly with the same technique with forthcoming 
HERA data, will be useful for the QCD analysis of present and future experiments at the Tevatron and the LHC. 
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